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^SJ Abstract. Second-harmonic-generation (SHG) in centrosymmetric material such as Si(lll) is 

t-H usually understood either from the phenomenology theory or more recently using the simplified 

CN bond-hyperpolarizability model (SBHM) [G. D. Powell, J. F. Wang, and D. E. Aspnes, Phys. 

Rev. B 65, 205320/1 (2002)]. Although SBHM is derived from a classical point of view, it 
^ has the advantage over the former that it gives especially for lower symmetry systems a clear 

• »"H physical picture and a more efficient explanation of how nonlinearity is generated. In this paper 

p\ we provide a step-by-step description of the SBHM in Si(lll) for the linear and second harmonic 

case. We present a SHG experiment of Si(lll) and show how it can be modelled by summing 

up the contribution of the fields produced by the anharmonic motion of the charges along the 

bonds. 
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1. Introduction 

Nonlinear optics is one of the most versatile tools to investigate surface properties, especially 
for a material class which posseses inversion symmetry. Serious theoretical efforts to explain 
nonlinearity in such systems has significantly progressed especially since the pioneering work 
of Bloembergen et.al. [1]. A nonlinear optical probe has several advantages over other surface 
probes because the material damage and contamination associated with charged particle beams 
is eliminated, all pressure ranges are accessible, insulators can be studied without the problem 
of charging effects, and buried interfaces are accessible owing to the large penetration depth of 
the optical radiation [2]. In centrosymmetric material such as silicon (Si) the bulk and surface 
nonlinear contributions are comparable [3] and changes in the surface properties such as surface 
deposition will significantly alter the nonlinear intensity profile [4]. 

Because centrosymmetric materials is important to present day technology (e.g. silicon 
semiconductor, silicon thin film sensors) understanding the surface property and the physical 
mechanism behind is critical, with much effort both experimental and theoretical having been 
invented. Perhaps, the most notable development in understanding the creation of an outgoing 
2uj wave or second-harmonic-generation (SHG) from centrosymmetric material surfaces and 
interfaces in the viewpoint of the phenomenological theory was performed by the group of Sipe 
et.al. [5] and further by Luepke et.al. [6]. The later group, calculated the nonlinear susceptibility 
for vicinal Si(lll) by fitting a fourier series to reproduce the four polarization SHG intensity 
(p-incoming p-outgoing, ps, sp, and ss) obtained from their experiments. Unfortunately, already 
for a single atom with a relatively low symmetry such as silicon, the phenomenological theory 
requires many -seemingly non independent- variables to be fitted thus blurring the physical 
insights and mechanism of SHG inside the material. 

As a response to this inefficiency, Aspnes et.al. [7] developed an ingenious method termed 
the simplified bond-hyperpolarizability model (SBHM) based on the somehow forgotten classical 
Ewald-Oseen extinction theory [8,9]. The theorem states that the electric field transmission and 
reflection can be obtained by direct superposition of dipoles inside the matherial rather than 
from macroscopic boundary calculation [10]. SBHM extends this classical view to cope with 
nonlinear optics by assuming that the SHG signal originates from the anharmonic motions of 
the charges along the bonds. Suprisingly, if the nonlinear field is far smaller than the linear field 
the calculation turns out to be much simpler because it does not lead to into an integral equation 
hence requiring self consistency check [11]. In this case it is possible to differentiate between 
the incoming driving fundamental field (w) and the outgoing SHG field (2w). As a result, the 
steps to obtain the far field contribution from the anharmonic bond radiation can be performed 
independently. Moreover, experimental data, e.g.Luepke's Si(lll) experiment, can be very well 
fitted using SBHM by only a handful parameters rather than the cumbersome fitting of various 
fresnel parameters using phenomenology theory. 

2. Basic Theory 

An excellent discussion of how to calculate the far field intensity for centrosymmtric material 
such as Si(lll) and Si(001) by assuming charge motion along their bonds already exists for 
vicinal Si sytems [6,12] and the interested reader is therefore suggested to refer to these papers 
for further information. However, for the sake of clarity, we will -although only briefly- show 
how SHG can be understood in the framework of the SBHM. 

In classical electrodynamics one can obtain the far field expression as a result of dipole (or 
oscillating charge) radiation. The vector potential A(f,t), which satisfy the Lorenz gauge [13]: 
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can be expressed in the form: 
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Assuming that the localized current j(r) and charge density p(r) as well as the vector potential 
A(f,t) can be written in harmonic form we can write eq.(2) in the far field approximation as: 
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Using the continuity equation: 

V • j(r') = iojp(r) 
we have for the spatial dependence vector potential the expression: 
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where p is the dipole moment: 



p = J d 3 r'r'p(r) 

The expressions for the fields according to the defined vector potential are: 
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The far field can thus be calculated by inserting eq.(5) in eq.(8) and then inserting eq.(8) in 
eq.(7). With the help of the identity {fix p) x n = p — (hh) ■ p we have: 



(9) 



here / is the identity tensor. The next step will be finding an expression for the classical 
microscopic dipole per unit volume (polarization) p. In this work we will use eq. (9) to calculate 
the reflected intensity of the linear and second-harmonic-generation. 



3. SBHM for Si(lll) 

In the Maxwell classical approach all the microscopic details in the inhomogenous nonlinear wave 
equation are included and hidden in the electric susceptibility. In the classical optic approach 
this susceptibilty can be estimated by constructing a Lorentz oscillator model with damping 
to obtain the nonlinear susceptibility terms [14]. It is therefore not suprising -but ingenious- 
to use a similar approach microscopically, by starting from a force equation to find expressions 
for the polarization. Furthermore it has been shown that even with the simplification that the 
nonlinear polarization source only occur along the bond direction, the second-harmonic far field 
from the anharmonic charge radiation is in very good agreement with experimental result for 
vicinal Si(lll). Here, we briefly follow the approach given by Ref. [7]. 



For a single Si(lll) atom there are four bonds which we denote by 61,62,63 and 64. In the 
bulk, the orientation of the bonds can be expressed in cartesian coordinates as: 
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with f3 = 109.47° is the angle between two bonds. Fig.l depicts a graphical image of the bond 
direction and the definition of the coordinate. Note that the four Si(lll) bonds are drawn as 
the red lines. 

Although the orientation is valid for the bulk, the even nonlinear harmonic dipole contribution 
is generally considered small because it is forbidden by inversion symmetry or parity argument. 
However we believe that the incoming field decay effects due to absorption may break the 
symmetry and produce dipolar effects even inside the bulk (spatial dispersion). At the interface, 
a dipole contribution is allowed because of symmetry breaking, but the orientation of the bonds 
may not be arranged as homogenously as in the bulk. Nonetheless, we follow the assumption 
in Ref. [7] that, based on statistical averaging the majority of the bond vector at the interface 
can be treated similar as within the bulk, thus the bond model of the bulk also holds for the 
interface. 

To model the reflected intensity vs rotation along the azimuthal direction (rotation in the 




Figure 1. Sketch of Si lattice ((111) surface) and definition of the bond direction in cartesian 
coordinates. 



x-y plane) we allow the bonds as a function of rotation via: 




bj((p) = bj ■ R<f, = bj ■ sin0 cos<^> (11) 



where the subscript j runs from 1 to 4. 

Having defined the bond vectors, we now analyze the motion of the charges. The force 
equation which describes this anharmonic motion along the bonds can be written in the form 
[7] taking the equilibrium position as zero: 

F = qjEbje~ tujt - K\x - k 2 x 2 - bx = mi (12) 

where q, m, x are the electron charge, mass, and its displacement from equilibrium, respectively 
and K\ and k 2 are the harmonic and anharmonic spring constants, and the term — bx is the 
common frictional loss in oscillation. Solving for Ax\ and Ax 2 by using the assumption that x 
can be written as x = xq + Axie~ lult Ax 2 e~ l2u)t gives for the lowest order of approximation: 




AX2 = A 2 -un ^—T- ( 13 ) 

therefore the linear polarization produced by each bond bj is 

Pij = Qj Ax i = «ij (h ■ E^j (14) 
whereas we have for the SHG the nonlinear polarization: 

p 2 j = qjAx 2 = a 2 j (bj ■ E^j (15) 

where a\ and a 2 are the microscopic polarizability and second order hyperpolarizability given 
by: 
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In Si(lll) surface, a.j=i is denoted as the "up" polarizability/hyperpolarizability a u whereas 
the three remaining have the same value due to symmetry and are called the "down" 
polarizability/hyperpolarizability ad- The total linear and second harmonic polarization 
produced by all the four bonds in Si(lll) considering azimuthal rotation are thus : 
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^ (P21 +P22 +P23 +P2A) = ^^2 a 2j (&j(0) " E) (18) 

In the following chapters, we will show how the SBHM can reproduce experimental results 
of second-harmonic-generation in Si(lll). 



4. SHG Experiment of Si(lll) 

The setup for the azimuthal SHG measurements is shown in Fig. 2. As a source of a radiation 
a compact femtosecond fiber laser system, working at a wavelength of 1560 nm and producing 
pulses with duration of 86 fs and with a repetition rate of 80 MHz, was used. After a variable 
attenuator for adjusting the power, a half-wave plate (HWP) determines the state of polarisation 
for the incident beam, which is then filtered by a band-pass filter (F), transmitting only the 
fundamental wavelength. The incidence radiation is focused onto the sample by the objective 
lens LI with a minimal focal waist down to 10 /im. 

Variable 
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Figure 2. Schematic sketch of the experimental setup for SHG measurements 
(Abbreviations: HWP: half-wave plate, F:filter, L:lens, G-T:Glan- Taylor polarizing prism, 
PMT:photomultiplier) . 

The incidence angle for the radiation was set to 45°. The sample itself is mounted on a 
motorized rotation stage. SH and fundamental radiation from the sample are collected by the 
lens L2 and are directed to a rotatable Glan- Taylor calcite polarizer that selects the desired 
polarisation direction. A Pellin Broca prism and a slit are used for spatially separating the 
SH radiation, which is focused onto a cooled photomultiplier tube (Hamamatsu R911P). For 
the acquisition of the SHG signal photon counting was used, with a dark count rate of only 30 
photons per minute. 

As a sample, a Si(lll) wafer with native oxide was used. In order to obtain a detectable 
SHG signal and not to damage the surface of the sample, the beam was slightly defocused, with 
the fiber laser operating at full power (350 mW). 

5. Si(lll) SHG experimental results and SBHM simulations 

Second harmonic generation from a Si(lll) face was obtained experimentally by plotting the 
intensity for various azimuthal angles (Si(lll) rotated by 360° in the x-y plane). All four 
intensity profiles (pp,ps, sp, ss) are aquired with the results depicted in Fig. 3(b). To obtain 



the total field using SBHM we insert expressions for the total polarization in eq. (17) and eq. 
(18) into the far field formula eq. (9), with n = k, where k is the outgoing wave vector (the unit 
vector in the direction of the observer) and is given by k = — cosOqx + sin 9qz. The linear and 
second harmonic far field can now be written as: 
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The term outside the paranthesis includes only constants and because SHG experiments are 
usually performed in arbitrary units the important terms are those given inside the paranthesis. 
Here g is the symbol for the polarization of the incoming wave and can be either p or s. 




Figure 3. Si(lll) linear and SHG intensity profiles vs azimuthal rotation in the x-y plane for 
various polarizations in arbitrary units of intensity, (a) The total linear pp-intensity yields 
a constant value independent of the azimutal rotation, (b) SHG intensity obtained from 
experiment and (c) SHG intensity simulation results using SBHM with a u =a,i=\. 



We next consider first eq. (17) and evaluate the intensity for four polarization direction: 
pp,ps, sp, and ss. The first letter in pp denotes the incoming and the second the outgoing wave. 
The far field from all four bonds in Si(lll) for each polarization can be obtained using eq.(10). 
For the linear case we have: 

Ew,pp = 7, (cos (Bo) x + sin (9 d ) z) [(-3a d cos (0j) cos (9 Q ) sin 2 (/3) + 3a d cos (2/3) sin (0;) sin (8 ))} 
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An interesting result from this calculation is that the total linear fields for all polarization 
is independent of the azimuthal rotation (ft, eventhough the field produced by each of the bond 
may depend on (ft. This can be seen in Fig. 3(a) where we take the pp-case as an example. 
Similar results are obtained for the ss-polarization. Therefore, orientation of the bonds due to 
rotation in the x-y plane does not effect the reflected linear intensity. 

More interesting results are obtained if we consider the second harmonic case. The total 
nonlinear polarization for all polarization modes are -in contrary to the linear case- a function 
of the azimuthal rotation of the crystal (ft. Fig. 3(b) shows our experimental result. The pp,ps, 
and ss-polarization (red line) showsa 6-fold dependence, with peaks at every 60° interval. The 
ps and ss-polarization have the same azimuthal behaviour with a peak shift of 30° for the 
pp-polarization. The ss-polarization, notably has a high (DC) offset value and broad peak full- 
wave- half-maximum (FWHM), but peaks still being detected. Also the height of the pp and 
sp-polarization peaks are quite similar with the sp polarization having a 3-fold symmetry. 

Remarkably, this seemingly complicated experimental result can be well explained using the 
SBHM, more explicitly by applying eq. (19) with the result depicted in Fig. 3(c). The formulas 
for the 4-polarization SHG far field from a Si(lll) structure can be seen in Ref. [7]: 
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with a slightly different result obtained for the pp-polarization: 

E20J, pp = (a u + 3ad cos 3 j3) sin 2 9i sin B — 3a d cos /3 cos 9i cos 9 Q sin 2 /3 sin 8i 
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By setting an arbitrary DC offsets which might be related to experimental noise, the model 
showed good agreement with the experiment, as is evident from the reproducable symmetry 
pattern with correct azimuthal position. In the simulation we have -for the sake of simplicity- set 
the values for the hyperpolarizabilities equal to unity and use for the SiO-Si interface incoming 
and outgoing angle of 29.5° (via snellius). The 3-fold sp symmetry (green line) is perfectly 
matching the experimental result as well as the ps-polarization blue line) with the correct relative 
peak difference between them (30°). 

The intensity profiles can be easily explained as an effect of anharmonic charge radiation along 
the bond direction. One can easily check this by calculating the field produced by the individual 
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bond contribution. For example, the 6-fold symmetry from the ss-polarization can be explained 
by each down bond contributing twice to the s-fundamental driving field when rotated by 360°. 
However, it has to be stated that for the pp-polarization SBHM gives a 3-fold symmetry and is 
in contradiction with the 6- fold symmetry given by experiment. We attribute this difference due 
to possible bulk dipole contribution because inside the bulk the incident angle is closer to the 
normal, and here SBHM also predicts a sixfold symmetry. This bulk dipole is still controversial 
and requires symmetry breaking in the form of a decaying field (absorption). Nevertheless, 
our simulation has shown that SHG measurement with SBHM can be used to investigate bond 
orientation at the surface of centrosymmetric materials in a more simple fashion without going 
to the complicated susceptibility tensor analysis (e.g. group theory). 

6. Summary 

Using the simplified bond hyperpolarizability model we show that the reflected linear intensity for 
all 4-polarization modes from a Si(lll) structure is independent of the azimuthal rotation on the 
x-y plane. For the second harmonic intensity, SBHM gives a good agreement with experiment, 
sucessfully predicting the correct azimuthal profiles and symmetries of the surfaces except for 
the pp-polarization which we believe requires SBHM to be extended to cope with bulk effects. 
This shows that SBHM can be used as a simple model in predicting the bond orientation of a 
Si(lll) wafer at the surface/interface and reconfirms second harmonic generation measurements 
as a sensitive non destructive method to investigate surface structures in material with inversion 
symmetry. 
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